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ABSTRACT. A deep far-infrared (12-lOOpm) survey is presented using over one 
thousand scans made of a 4-6 deg 2 field at the north ecliptic pole by the Infra- 
red Astronomical Satellite (IRAS). Point sources from this survey are up to 
one hundred times fainter than the IRAS point source catalog at 12 and 25pm, 
and up to ten times fainter at 60 and lOOym. The 12 and 25ym maps are instru- 
mental noise-limited, and the 60 and 100pm maps are confusion noise-limited. 

The majority of the 12ym point sources are stars within the Milky Way. The 
25ym sources are composed almost equally of stars and galaxies. About 80 
percent of the 60ym sources correspond to galaxies on Palomar Observatory Sky 
Survey (POSS) enlargements. The remaining 20 percent are probably galaxies 
below the POSS detection limit. The differential source counts are presented 
and compared with theoretical models. The 12ym source counts are found to be 
consistent with what is predicted by the Bahcall and Soneira Standard Galaxy 
Model (1980) using the B-V-12ym colors of stars without ci rcumstel lar dust 
shells given by Waters, Cote and Aumann (1986). The 60ym source counts are 
inconsistent with those predicted for a uniformly distributed, nonevolving 
universe. The implications are briefly discussed. A detailed description of 
the survey appears in Hacking and Houck (1986). A discussion of the scientific 
implications of the source counts and colors are discussed in Hacking, Condon, 
and Houck (1986) . 


INTRODUCTION 

There are two reasons for conducting an infrared survey to lower flux 
densities. The first is to detect trends in sources observed at higher flux 
densities. The changes might be detected in the shape of the luminosity 
function (i.e., changes in relative distribution of sources vs. luminosity) or 
the position of the luminosity function in the density — luminosity plane 
(e.g., changes in the space density of sources or luminosity evolution at 
cosmological redshifts). The second is to search for new classes or types of 
objects. An obvious example would be the detection of a nearby , isolated brown 
dwarf. It is also possible that a new type of source, as yet unthought of, 
might be detected. 


THE DATA 

The survey is located at the north ecliptic pole, in the constellation of 
Draco, at 18 h , +66°; Jill = 97°, b 11 = 30°, and covers 4.3 square degrees 
at 12 and 25pm and 6.3 square degrees at 60 and 100pm. 

There are two types of data used in the survey: 488 scans from the IRAS 
all-sky survey and 141 pointed observations. The survey scans cover the entire 
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field approximately uniformly. Each survey scan map is 1/2° x up to 4° and 
consists of data from a single scan of the satellite across a portion of the 
field. The pointed observations are centered on NGC 6543, which is located 
very near the ecliptic pole and was used as a calibrator for the IRAS mission. 
Each pointed observation consists of three to six scans made at half the IRAS 
survey scan rate and results in a map of dimensions 1/2° x 1 to 1.75°. Since 
NGC 6543 was the target of the pointed observations, only the central three 
square degrees is covered by them, and the coverage is nonuniform, with the 
region around NGC 6543 (near the center of the field) observed the most. Over 
one thousand scans were combined to form the final maps, making this the 
deepest infrared survey to date. At 12 and 25ym it is likely to remain so 
until the next generation infrared satellite can be flown (e.g., SIRTF or 
ISO). At 60 and lOOym, it may prove possible to probe slightly deeper in the 
field (and a few other smaller select fields) using different filtering 
techniques. Point sources as faint as 5-10 mJy are in the 12 and 25ym maps 
with SNR > 3-5 (a factor of 50 to 100 times fainter than the IRAS point source 
catalog), 50 mJy at 60ym (ten times fainter than the IRAS PSC), and 100-200 mJy 
at lOOym (five to ten times fainter than the IRAS PSC). 

Two images of each survey scan and each pointed observation were made: 

1) A ful 1 -intensity image was made that approximates the true appearance of 
the infrared sky; 2) a point-source filtered image was made that improves the 
si gnal -to-noi se ratio of point sources. Several preliminary maps were made 
from subsets of the data to check for moving or spurious sources. All of the 
images were then summed to form final intensity and point source filtered 
images. 

The intensity maps indicate that there is a great deal of extended emis- 
sion in the field at 60 and lOOym. The 60ym point source filtered map contains 
very little extended emission. The lOOym point source filtered map is contami- 
nated with extended emission, however. This is due to the larger angular size 
of the lOOym point spread function that is similar in size to much of the 
structure present in the extended emission. For this reason, lOOym fluxes are 
quoted only for sources also detected at 60ym, where contamination from the 
extended emission is less likely. 


RESULTS 

The point source density in this field at 12 and 25ym is quite low 
coup a red to the effective beam size (10-15 sources per square degree brighter 
than 20 mJy), so that the 12 and 25ym surveys are limited by instrumental 
noi se. 

At 60ym, the effective beam size and the source density (20 sources per 
square degree brighter than 50 mJy) are larger, resulting in a confusion 
noise-limited survey (a con f = 20 mJy). Extended emission and confusion 
prevented an unbiased lOOym survey. 

No moving or spurious (single event) sources were found in any of the 
preliminary maps. 

The point source filtered maps were used to extract the point source 
sample. The point sources were compared with the Smithsonian Ast rophysi cal 
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Observatory Star Catalog (1966) and the Palomar Observatory Sky Survey (POSS) 
plates. A summary of the comparison is presented in Table 1. 


TABLE 1 

- Summary 

of Point 

Sources from 

Deep Survey 


Total 

Stars 

Gal axi es 

Unidenti fied 

X 

Numbe r a 

% 

% 

% 

12pm 

46 

89 

11 

0 

25pm 

36 

50 

47 

3b 

60pm 

99 

0 

~80 

~20C 


a Excludi ng NGC 6643 
b Probably a distant Seyfert galaxy 
c Probably galaxies belcw POSS detection limit 


The distribution of sources vs. flux density was modeled, assuni ng a power 
law of the form: 


where n>p v is the number of sources with flux density greater than F v per 
square degree; A and a are constants. Values of a of U.76±0.11 and 1.79±0.17 
were found for the 12pm (mostly stars) and 60pm (ext ragal acti c sources) point 
source samples, using the method of maximum likelihood. The Kolmogorov test 
revealed that the power law distribution assumed for the source counts is 
consistent with the data at 12 and 60pm. Since the 25pm sample is composed of 
both stars and galaxies, the 25pm distribution is a linear combination of the 
12 and 60pm distributions. 


DISCUSSION 

The 12pm source counts are consistent with star counts predicted by the 
Bahcall and Soneira Standard Galaxy Model*, using the B-V-12pm colors given 
by Waters, Cote, and Aumann (1986) for stars without ci rcumstel lar dust 
shells. In particular, the total number (33) and the value of a, 2 (0.78) 

is successfully predicted (see Hacking, Condon, and Houck 1986 foraetails). 
Figure 1 shows differential star counts extrapolated to very low flux densities 
at the north galactic pole (not covered in this survey) as predicted by the 
model. The 12pm counts above ~10Jy are dominated by stars with ci rcumstel lar 
dust shells (Hacking, et al . 1985), which are not included in Figure 1. The 
expected galaxy counts at 12pm are also plotted in Figure 1. As can be seen, 
galaxies should dominate the 12pm point source counts below ~10 mJy at the 
galactic poles. 

The 60pm source counts were modeled in a similar fashion to the radio 

*See Hacking, Condon, and Houck (1986) for a description of what was used, 
and Bahcall and Soneira (1980), for a detailed description of the model. 
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Figure 1 - Predicted 12ym differential star counts at the north galactic pole 
using the Bahcall and Sonei ra Standard Galaxy Model (1980) and 
B-V-12ym star colors from Waters, Cote, and Aumann (1986). The log 
of the number of sources per square degree, per decade of flux 
density is plotted versus log flux density (in Janskys). The 
differential counts for all stars are given by the solid line. 

In addition, the differential counts for main sequence and giant 
stars are each shown for the disk and spheroidal components. The 
expected galaxy counts at 12ym are also shown. 

source counts using the model of Condon (1984) with and without evolution, 
using the 60ym luminosity function and color-luminosity data given by Soifer, 
et al . (1986). The predicted counts for an evolutionary model and a 
nonevolving model are compared with the observed counts in Figure 2. Although 
the source counts favor the evolving model, they do not constrain the type of 
evolution that may have been detected. In addition, evolution is not the only 
explanation for the discrepancy between the observed source counts and the 
source counts predicted for a nonevolving universe. In particular, a galaxy 
cluster in the direction of this survey with a redshift of ~0.2 could produce 
this effect. Clearly, redshift data are needed to resolve the nature of this 
di screpang . 
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Figure 2 - Predicted 60ym differential galaxy counts with and without 

evolution, using the luminosity function and luninosity-color data 
from Soifer et al . (1986). Tick marks i ndi cate medi an redshi fts of 
the differential counts. The data from the deep survey are plotted 
as filled circles, and the data from the IRAS Point Source Catalog 
with I b | > 50° are plotted as filled squares. Much of the increase 
in the Source counts above 10 Jy is due to the Virgo cluster. The 
luminosity function has been shifted down by 15% so that the 
non-evolutionary model agrees with the differential counts at 1 
Jy. See Hacking, Condon, and Houck (1986) for details. 
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DISCUSSION 
WYNN- WILLI AMS: 

Are there other IRAS fields that can be analyzed to this depth, or have we hit the limit of what is possible 
with IRAS? 

HACKING: 

This is probably as deep as we can go at 12 and 25/; m. There are smaller deep fields that are almost as 
deep as this field at 60 and 100/im. In addition, it may prove possible to reduce the confusion noise in this 
field using different filtering techniques, allowing a study up to a factor of 2 fainter in this field. 

WEEDMAN: 

Your counts at 60//rn agree precisely with expectations from the average luminosity function now avail- 
able. The very close agreement must be fortuitous, but it is pleasing that all of the data are coming together 
so consistently. 
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